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Statement  of  the  Problem  Solved 

Challenges  addressed  by  the  research  include:  (i)  improving  our  understanding  of  the  many- 
electron  quantum  mechanisms  by  which  nature  uses  strong  electron  correlation  for  efficient 
energy  transfer,  particularly  in  photosynthesis  and  bioluminescence,  (ii)  providing  an  innovative 
paradigm  for  energy  transfer  in  photovoltaic  materials  by  which  new  levels  of  solar  efficiency 
are  achieved  through  the  use  of  strong  electron  correlation  and  entanglement,  (iii)  enhancing 
two-electron  reduced-density-matrix  (2-RDM)-based  electronic-structure  methods  that 
significantly  expand  the  range  of  strongly  correlated  molecular  systems  that  can  be  studied  with 
applications  throughout  science  and  engineering,  and  (iv)  developing  a  new  generation  of  large- 
scale,  parallel  algorithms  for  performing  semidefinite  programming  with  applications  to 
problems  in  engineering,  computer  science,  statistics,  finance  and  economics. 

Summary  of  the  Most  Important  Results 

Significant  research  was  supported  under  a  prior  grant  from  the  Army  Research  Office. 
The  research  produced  33  published  papers1"33  and  partially  supported  the  training  of  12  graduate 
students,  5  undergraduate  students,  and  1  high  school  student.  The  high  school  student  Claire 
Liu  was  recently  named  one  of  the  semifinalists  in  the  2014  Intel  Science  Talent  Search  for  her 
research  project,  “The  Effect  of  Electron  Correlation  on  Efficient  Energy  Transfer  in  a  Synthetic 
Chromophore  System.”  Highlights  of  the  research  include: 

1)  Development,  Application,  and  Distribution  of  the  Parametric  2-RDM  Method: 
Under  the  ARO  grant  we  developed  the  parametric  2-RDM  method  (p2-RDM)  in  which  the  2- 
RDM  is  parameterized  by  single  and  double  excitations.  The  p2-RDM  method  has  a  smaller 
pre-factor  in  its  computational  scaling  than  coupled  cluster  singles-doubles  with  an  accuracy 
approaching  wave  function  methods  with  triple  excitations,  especially  in  the  presence  of  multi¬ 
reference  correlation.  Applications  have  been  made  to  strongly  correlated  transition  states  and 
intermediates  including  the  rotational  transition  state  between  the  cis  and  trans  diazene  as  well  as 
the  diradical  isomers  of  the  recently  discovered  olympicene  molecule.  The  method  has  also  been 
applied  to  showing  the  importance  of  multi-reference  correlation  on  small  energy  differences 
which  arise  in  determining  the  existence  and  stability  of  oxywater  and  the  relative  energies  and 
stability  of  the  cage  and  prism  isomers  of  water  hexamer.  The  p2-RDM  method  has  made 
available  for  download  from  the  research  group’s  website. 

2)  Prototype  Low -scaling  2-RDM  Method:  Using  tensor  factorizations  of  the  double¬ 
excitation  and  Hamiltonian  matrices,  we  developed  a  prototype,  low-scaling  2-RDM  method  that 
scales  computationally  as  0(r4)  where  r  is  the  number  of  orbitals  in  the  one  electron  basis  set. 
Further  development  of  the  method  has  the  potential  to  provide  an  ab  initio  alternative  to  density 
functional  theory. 

3)  Complete  Set  of  N-representability  Conditions:  The  2-RDM  must  be  constrained  for 
its  two  electrons  to  represent  an  A-electron  quantum  system.  These  constraints  are  known  as  N- 
representability  conditions.  With  ARO  support  we  presented  a  constructive  solution  to  the  N- 
representability  problem:  a  full  characterization  of  the  conditions  for  constraining  the  two- 
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electron  reduced  density  matrix  to  represent  an  iV-electron  density  matrix.  Previously  known 
conditions,  while  rigorous,  were  incomplete.  We  derived  a  hierarchy  of  constraints  built  upon  (i) 
the  bipolar  theorem  and  (ii)  tensor  decompositions  of  model  Hamiltonians.  Existing  conditions 
(D,  Q,  G,  Tl,  and  T2),  known  classical  conditions,  and  new  conditions  appear  naturally.  This 
theoretical  advance  gives  completeness  to  the  foundation  of  2-RDM  quantum  theory.  Conditions 
beyond  T2  will  eventually  allow  for  the  computation  of  molecular  quantum  systems  in 
polynomial  time  with  unprecedented  accuracy. 

4)  Efficient  Algorithm  for  Large-scale  Semidefinite  Programming:  We  developed  a 
more  efficient  approach  to  adding  Wrepresentability  conditions  to  2-RDM  calculations  based  on 
a  boundary-point  algorithm  for  semidefinite  programming.  The  current  parallel  version  of  the 
boundary-point  algorithm  is  100  to  1000  times  faster  than  a  previous  first-order  algorithm  for 
semidefinite  programming  developed  by  our  group  for  electronic  structure  calculations.  The 
algorithm  has  been  adopted  by  others  for  2-RDM  calculations  as  well  as  other  applications  in 
science  and  engineering. 

5)  Functional  Subsystems  in  Photosynthetic  Light  Harvesting:  The  Fenna-Matthews- 
Olson  (FMO)  antennae  complex,  responsible  for  light  harvesting  in  green  sulfur  bacteria, 
consists  of  three  monomers,  each  with  seven  chromophores.  Skochdopole  and  Mazziotti  showed 
that  multiple  subsystems  of  the  seven  chromophores  can  transfer  energy  from  either 
chromophore  1  or  6  to  the  reaction  center  with  an  efficiency  matching  or  in  many  cases 
exceeding  that  of  the  full  seven  chromophore  system.  In  the  FMO  complex,  these  functional 
subsystems  support  multiple  quantum  pathways  for  efficient  energy  transfer  that  provide  a  built- 
in  quantum  redundancy.  There  are  many  instances  of  redundancy  in  nature,  providing  reliability 
and  protection,  and  in  photosynthetic  light  harvesting  this  quantum  redundancy  provides 
protection  against  the  temporary  or  permanent  loss  of  one  or  more  chromophores.  The  complete 
characterization  of  functional  subsystems  within  the  FMO  complex  offers  a  detailed  map  of  the 
energy  flow  within  the  FMO  complex,  which  has  potential  applications  to  the  design  of  more 
efficient  photovoltaic  devices. 

6)  Strong  Correlation  in  Light  Harvesting:  We  generalized  the  single-electron/exciton 
models  to  a  multi-electron/exciton  model  that  explicitly  shows  the  effects  of  enhanced  electron 
correlation  within  chromophores  on  the  efficiency  of  energy  transfer.  Importantly,  within  the 
model  the  addition  of  electron  correlation  within  the  chromophore  enhances  energy  transfer  by 
as  much  as  100%. 

7)  Measurement-driven  Reconstruction  of  Many-particle  Quantum  Processes: 
Quantum  measurements  provide  a  trove  of  information  about  a  quantum  system  or  process 
without  solution  of  the  Schrodinger  equation,  and  in  principle,  the  associated  density  matrix  is  a 
function  of  these  measurements.  Inversion  of  the  measurements  can  produce  an  estimate  of  the 
density  matrix,  but  this  estimate  may  be  unphysical,  especially  when  the  measurements  are  noisy 
or  incomplete.  We  have  recently  developed  a  general  approach  based  on  semidefinite 
programming  for  reconstructing  the  density  matrix  from  quantum  measurements  which  leads 
naturally  to  nonnegative  solutions,  a  critical  attribute  of  physically  realistic  solutions.  We  have 
applied  the  methodology  to  reconstructing  the  time -dependent  quantum  process  of  exciton 
transfer  in  a  photosynthetic  light-harvesting  complex.  The  research  also  has  broad  applications  in 
advanced  sensing  and  image  reconstruction. 

8)  Metrics  for  Assessing  Electron  Correlation  and  Entanglement:  Reliable  but 
computable  metrics  for  assessing  the  degree  of  electron  correlation  and  entanglement  are  critical 
for  studying  strongly  correlated  molecules  and  materials.  Reduced  density  matrices  (RDMs)  are 
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natural  variables  for  measuring  the  degree  of  correlation  between  electrons.  We  studied  the  use 
of  the  2-RDM  and  parts  of  the  2-RDM  for  the  measurement  of  electron  correlation  and 
entanglement  in  the  photosynthetic  light  harvesting  of  green-sulfur  bacteria.  We  found  that  the 
squared  Frobenius  norm  of  the  cumulant  2-RDM  very  accurately  represents  the  entanglement  in 
this  complicated  time-dependent  process. 

8)  Strong  Correlation  in  Molecular  Systems  and  Processes:  The  2-RDM  methods  were 
applied  to  computing  the  energies  and  properties  of  strongly  correlated  molecular  systems  and 
processes  including  chains  (acenes)  and  sheets  of  polyaromatic  hydrocarbons  and  the  reaction  of 
firefly  luciferin  in  bioluminescence.  Calculations  of  both  ground-  and  excited-state  energies 
were  performed  at  the  conical  intersections  in  the  chemiluminescence  of  dioxetanone  and  the 
photoexcited  tautomerization  of  vinyl  alcohol  to  acetylaldehyde. 
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